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1 3.  ABSTRACT  ( Maximum  200  words) 

We  present  numerical  simulations  of  the  collision  and  subsequent  interaction  of  two  initially  orthogonal,  twisted,  force-free  field 
magnetic  fluxtubes.  The  simulations  were  carried  out  using  a  new  three-dimensional  explicit  parallelized  Fourier  collocation 
algorithm  for  solving  the  viscoresistive  equations  of  compressible  magnetohydrodynamics.  It  is  found  that,  under  a  wide  range  of 
conditions,  the  fluxtubes  can  "tunnel"  through  each  other.  Two  key  conditions  must  be  satisfied  for  tunneling  to  occur:  the  magnetic 
field  must  be  highly  twisted  with  a  field  line  pitch  >  >  1,  and  the  magnetic  Lundquist  number  must  be  somewhat  large,  ^  2880.  this 
tunneling  behavior  has  not  been  seen  previously  in  studies  of  either  vortex  tube  or  magnetic  fluxtube  interactions.  An  examination  of 
magnetic  field  lines  shows  that  tunneling  is  due  to  a  double-reconnection  mechanism.  Initially  orthogonal  field  lines  reconnect  at  two 
specific  locations,  exchange  interacting  sections  and  "pass"  through  each  other.  The  implications  of  these  results  for  solar  and  space 
plasmas  are  discussed. 
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MAGNETIC  FLUXTUBE  TUNNELING 

Introduction 

Models  of  magnetic  reconnection  in  the  Sun’s  interior  and  atmosphere  usually  begin 
with  a  purely  two-dimensional  geometry.  However,  the  magnetic  field  at  the  solar  pho¬ 
tosphere  is  observed  to  be  organized  into  isolated  flux  bundles  1,  a  structure  which  must 
continue  into  the  corona.  In  addition,  photospheric  motions  “wind  up”  this  magnetic  field, 
producing  twisted  fluxtubes  2  3.  Interaction  between  fluxtubes  has  been  proposed  as  a 
mechanism  for  coronal  heating  4  5  6  and  might  be  the  origin  of  the  fine-scale  temporal 
variability  of  hard  X-ray  and  microwave  emission  observed  in  two- ribbon  flaxes  7  8.  Flux¬ 
tubes  axe  widely  believed  to  be  the  dominant  magnetic  structure  in  the  convection  zone, 

9  Coronal  mass  ejections  have  been  identified  in  the  interplanetary  medium  as  flux  ropes, 

10  which  are  thought  to  be  essential  ingredients  for  reconnection  at  the  magnetopause  11 
and  magnetotail  12 .  Therefore,  a  key  issue  for  understanding  many  important  phenomena 
in  solar  and  space  physics  is  the  nature  of  fluxtube  interaction. 

We  are  using  numerical  simulations  to  investigate  the  basic  physics  of  magnetic  flux¬ 
tube  collision  and  reconnection.  Our  explicit,  Fourier  collocation  algorithm,  which  is  de¬ 
scribed  in  detail  elsewhere  13  14 ,  solves  the  three-dimensional,  compressible,  dissipative 
magnetofluid  equations  in  a  dimensionless  form14.  The  geometry  used  is  that  of  a  triply 
periodic  cube  with  sides  equal  to  2i r,  making  a  Fourier  spectral  method  the  optimal  choice 
for  spatial  discretization.  The  results  described  here  were  computed  by  a  parallelized  ver¬ 
sion  of  our  code  recently  implemented  on  the  256  processor  NRL  TMC  CM5E  15 .  A  typical 
resolution  for  the  runs  in  this  paper  is  1283  Fourier  modes,  requiring  approximately  5  s 
per  timestep  and  1.05  GB  of  parallel  memory. 
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The  initial  conditions  consist  of  two  orthogonal  fluxtubes  and  a  flow  field  which  drives 
them  together.  Each  of  the  tubes  is  initialized  using  the  Gold-Hoyle  model  of  a  uniformly 
twisted,  cylindrical,  force-free  magnetic  field  (B  =  (Bx,By,Bz))  16  ,  viz., 


B0br  sin  <f>  n  Bobrcos<f>  n  B0 

1  +  62r2  ’  y  =  1  +  62r2  5  Bz  =  1  +  6V2  ’ 


(1) 


where  Bq  —  4,  and  r  and  <j>  are  the  radial  and  cylindrical  coordinates  of  the  fluxtube. 
The  parameter  b  measures  the  field  line  pitch  z.e.,  b  =  d<f>/dz.  In  order  to  maintain  ideal 
equilibrium,  the  uniform  gas  pressure  (p)  outside  the  tube  is  set  top  =  p0+2Bq/(1  +  62  B2), 
where  p0  =  20/3  is  the  gas  pressure  inside  the  tube,  and  R  =  ll7r/48  is  the  fluxtube  radius. 
To  ameliorate  the  Gibbs  phenomenon  due  to  the  Fourier  series  discretization  of  the  sharp 
cutoff  at  r  —  R,  we  pass  these  initial  conditions  through  a  raised  cosine  filter  17 .  The 
initial  density  is  uniform,  p  =  1.  These  values  for  B0  and  po  yield  a  plasma  /3  =  0.42  at  the 
fluxtube  axis.  We  center  one  tube  at  ( x  =  57r/4,  z  =  7r).  A  second  tube,  rotated  by  7r/2,  is 
centered  at  (x  =  37r/4,  y  =  37r/4)  to  break  the  symmetry.  There  axe  four  physically  distinct 
relative  orientations  for  the  fluxtubes,  depending  on  the  choice  of  the  axial  and  azimuthal 
magnetic  field  in  each  of  the  tubes.  In  all  our  simulations  the  twist  and  orientation  are 
chosen  so  as  to  maximize  the  possibility  of  reconnection  —  where  they  collide,  the  axial 
(azimuthal)  field  of  the  horizontal  tube  in  the  collision  region  is  directed  opposite  to  the 
azimuthal  (axial)  field  of  the  vertical  tube,  as  shown  in  Figure  1. 

The  two  fluxtubes  are  driven  together  by  an  initial  velocity  field  given  by: 


v(x,  y,  2,  t  =  0)  =  ev  —  sin  x(cos  y  +  cos  z)ex  +  cos  x(sin  ye  y  +  sin  zez ) 


(2) 
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The  velocity  amplitude  is  chosen  to  be  2.5%  of  the  Alfven  speed  at  tube  center,  ev  =  0.1, 
which  is  equivalent  to  4.25%  of  the  sound  speed  outside  the  tubes.  It  should  be  emphasized 
that  we  do  not  impose  any  subsequent  driver  on  the  system,  so  that  the  fluxtubes  evolve 
freely.  Uniform  and  isotropic  viscosity  and  resistivity  axe  used,  and  are  chosen  so  that  the 
resistive  and  viscous  Lundquist  numbers  are  equal.  The  two  parameters  that  are  varied 
in  our  system  are  the  field  line  pitch  or  twist,  6 ,  and  the  Lundquist  number,  S .  Low  and 
high  twist  (6  =  1  and  10),  and  low  and  high  Lundquist  number  ( S  —  576  and  2880)  runs 
have  been  performed. 


Results 

In  all  cases  we  find  that  the  fluxtubes  collide  and  begin  to  reconnect  as  a  result  of 
their  initial  velocity.  The  standard  picture  for  this  process  is  that  the  magnetic  field  of 
one  tube  opposes  the  motion  of  the  other  tube.  A  thin  current  layer  builds  up  between 
the  two  tubes,  and  if  the  magnetofluid  is  resistive  and  the  orientation  of  the  field  fines 
favorable,  reconnection  occurs  at  the  Petschek  rate  18 ,  the  Sweet-Parker  rate  9,  or  the 
linear  tearing  rate  7 .  The  reconnected  field  fines  will  then  pull  away  from  the  reconnection 
site:  this  is  the  basis  of  particle  acceleration  by  the  “slingshot  effect”  19 .  The  expectation 
is  that  topologically  the  magnetic  fluxtubes  will  reconnect  like  orthogonal  vortex  tubes  20 
so  that  two  initially  straight  orthogonal  tubes  exchange  halves  to  form  two  tubes  bent  at 
right-angles. 

In  contrast  to  this  picture,  we  find  that  for  low  twist  reconnection  between  the  fluxtubes 
is  not  substantial,  even  for  low  Lundquist  numbers.  The  evolution  is  primarily  an  elastic 
collision  of  the  fluxtubes  with  little  energy  released  to  the  plasma.  This  result  is  important 
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to  models  of  solar  activity,  because  coronal  fluxtubes  are  believed  to  have  low  twist.  Our 
results  demonstrate  that,  at  least  for  modest  collision  velocities,  fluxtube  reconnection  is 
much  more  difficult  to  accomplish  than  generally  is  assumed. 

We  then  attempted  to  obtain  more  reconnection  by  increasing  the  twist  of  the  fluxtubes 
to  6  =  10  and  by  using  a  low  Lundquist  number  ( S  =  576).  Two  effects  can  enhance 
reconnection  for  higher  twist.  First,  two  null  points,  at  which  opposing  field  lines  are 
exactly  anti-parallel,  can  form  in  the  interaction  region.  It  is  widely  believed  that  rapid 
reconnection  is  favored  at  such  nulls  21 .  Second,  as  we  show  below,  the  reconnected  field 
lines  are  highly  tangled  and  tend  to  pull  the  tubes  together,  thereby  enhancing  their 
interaction.  This  tangling  appears  to  be  the  dominant  effect.  We  find  that  the  high-6,  low¬ 
s'  case  does,  indeed,  behave  like  the  standard  vortex-tube  reconnection  picture  described 
above  15 . 

Since  low  Lundquist  numbers  generally  are  not  relevant  to  space  plasmas,  we  then 
considered  the  effect  of  increasing  S;  our  expectation  being  that,  for  sufficiently  high  S, 
the  reconnection  rate  would  decrease  to  the  point  that  the  tubes  behave  like  the  low  twist 
case  and  simply  bounce  off  each  other.  The  actual  evolution  was  a  complete  surprise  —  for 
S  =  2880  the  tubes  apparently  pass  right  through  each  other.  This  result  is  demonstrated 
in  Figure  2,  which  shows  isosurfaces  of  the  magnetic  field  magnitude  at  six  times  during 
the  run.  The  isosurfaces  are  taken  at  half  the  maximum  field  magnitude  value  at  each 
time.  Initially,  the  horizontal  fluxtube  is  behind  the  vertical  fluxtube.  The  effect  of  the 
first  reconnections  becomes  clear  as  the  reconnected  field  lines  exert  a  torque  on  the  tubes, 
giving  them  a  kinked  appearance.  Later,  the  singly  connected  field  lines  dominate,  as  the 
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tubes  begin  to  exchange  halves.  Then  the  doubly  reconnected  field  dominates,  with  the 
central  field  lines  pulling  along  the  rest  of  the  field  lines.  At  the  final  time,  the  vertical 
tube  is  now  in  front  of  the  horizontal  tube,  in  contrast  with  the  initial  state. 

The  origin  of  tunneling  can  be  clarified  by  close  inspection  of  the  magnetic  connectivity. 
Figure  1  shows  one  winding  of  two  representative  field  lines  before  reconnection.  Figure 
3  shows  the  effect  of  the  first  reconnection,  (for  clarity  only  one  of  the  reconnected  field 
lines  is  shown).  Since  there  are  two  reconnection  regions,  each  field  line  tends  to  reconnect 
twice.  The  lines  exchange  halves  to  form  two  right-angle  lines,  which  cannot  simply  pull 
away  from  the  reconnection  site  as  in  the  standard  picture,  however,  because  they  are 
wrapped  around  both  fluxtubes.  In  fact,  this  tangling  tends  to  pull  the  tubes  together, 
producing  more  reconnection.  Figure  4  shows  the  effects  of  the  second  reconnection.  The 
vertical  field  line  passes  around  the  horizontal  tube,  while  a  similar  horizontal  field  line 
passes  around  the  vertical  tube.  Such  field  fines  are  now  free  to  continue  on  their  initial 
trajectory  without  further  interaction.  This  is  the  basic  mechanism  of  magnetic  tunneling; 
two  reconnections  at  approximately  fixed,  diagonally  opposite  points  in  the  collision  region, 
(see  Figure  5),  allow  vertical  and  horizontal  field  fines  to  exchange  central  portions  and, 
thereby,  pass  through  each  other. 

Discussion 

Several  interesting  features  and  unanswered  questions  have  been  raised  by  our  results. 
First,  it  is  evident  from  Figure  2  that  the  post-reconnection  fluxtubes  have  very  com¬ 
plex  internal  structure.  In  the  initial  Gold-Hoyle  model,  all  field  fines  fie  on  well-defined 
cylindrical  flux  surfaces.  Some  of  the  initial  twist  is  lost,  however,  during  the  double  re- 
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connection  process  —  at  least  one  winding.  For  high  S,  magnetic  helicity  is  expected  to 
be  approximately  conserved  22 ;  hence,  the  helicity  in  the  twist  is  transferred  to  internal 
wrappings  of  the  field  lines  about  each  other.  The  conclusion  is  that,  as  in  the  case  of 
2.5D  studies  23 ,  reconnection  tends  to  produce  a  chaotic  final  field  topology.  These  results 
imply  that  magnetic  fluxtubes  which  erupt  through  the  solar  surface  also  axe  likely  to  have 
a  complex  fine-scale  internal  structure. 

One  fascinating  aspect  of  tunneling  is  the  process  by  which  the  central  field  lines  near 
the  tube  axis  pass  through  each  other.  The  double  reconnection  mechanism  described 
above  should  be  inapplicable  to  them  since  these  field  lines  have  vanishing  twist.  In  fact, 
they  also  undergo  a  double  reconnection  similar  to  the  highly  twisted  lines.  Figure  2  shows 
that,  during  the  height  of  the  reconnection  phase,  both  fluxtubes  are  deformed  so  that  they 
run  diagonally  through  the  reconnection  sites.  We  believe  that  this  deformation  is  due  to 
field  lines,  such  as  those  in  Figure  4,  which  have  already  tunneled  and  exert  a  torque  on 
the  remaining  field,  causing  the  central  field  lines  to  become  parallel  in  the  interaction 
region.  An  additional  force  pulls  the  axis  field  lines  out  of  their  respective  vertical  planes 
and  toward  the  two  reconnection  sites  where  they  also  reconnect  and  exchange  central 
portions,  as  seen  in  Figure  5.  The  origin  of  this  force  is  unclear  to  us,  at  present.  We  also 
do  not  know  whether  the  observed  evolution  reflects  some  general  physical  property  of 
magnetic  fluxtubes  or  is  simply  a  peculiarity  of  the  system  under  consideration.  It  should 
be  noted  that,  for  the  high  twist  case,  the  energy  in  the  azimuthal  field,  27r  f  B^2  rdr , 
is  over  3  times  the  axial  field  energy.  This  is  another  reason  for  expecting  tunneling  to 
require  a  large  twist.  Unless  the  twist  component  dominates,  it  cannot  deform  the  axial 
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component  sufficiently  for  that  component  to  also  undergo  double  reconnection.  A  further 
possibly  important  point  is  that  high  twist  fluxtubes  are  expected  to  be  kink  unstable. 
Since  the  Gold-Hoyle  model  is  metastable  to  first  order,  the  kink  mode  growth  rate  is 
small,  24but  it  may  be  contributing  to  the  deformation  of  the  axial  field. 

In  the  context  of  solar  and  space  plasmas,  the  key  question  is  the  range  of  applicability 
of  our  results  to  the  very  large  Lundquist  number  regime.  We  have  two  main  conclusions 
-  for  modest  collision  velocities,  low  twist  orthogonal  tubes  bounce  and  high  twist  tubes 
tunnel.  The  first  result  should  be  even  more  valid  at  higher  S;  hence,  fluxtube  reconnection 
is  unlikely  to  be  common  in  the  solar  corona  unless  the  tubes  are  driven  together,  or  are 
nearly  anti-parallel.  We  believe  that  the  second  result  also  holds  for  higher  S  because 
the  rate  of  reconnection  depends  very  weakly  on  5,  but  this  is  only  a  conjecture.  We  are 
currently  exploring  these  issues. 
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FIGURE  2.  Isosurfaces  of  |B|  at  several  times  during  the  run.  The  upper  left  corner 
shows  the  result  for  t  —  7.8.  To  the  right  are  the  results  for  t  =  21.1  and  30  6  The 
iower  left  corner  shows  the  result  for  t  =  40.8.  To  its  right  are  the  results  for  t  =  52  9 
and  64.2.  The  isosurfaces  are  chosen  to  equal  ||B|max  at  each  time. 
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